arXiv:1506.05720vl [nlin.CD] 18Jun2015 


Multifractality of quantum wave functions in the presence of perturbations 


R. Dubertrand, 1 ’ 2 ' 3 I. Garcla-Mata, 4 - 5 B. Georgeot, 1,2 O. Giraud, 6 G. Lemarie, 1 ’ 2 and J. Martin 3 

'Universite de Toulouse; UPS; Laboratoire de Physique Theorique (IRSAMC); F-31062 Toulouse, France 
2 CNRS; LPT (IRSAMC); UMR 5152, F-31062 Toulouse, France 
3 Institut de Physique Nucleaire, Atomique et de Spectroscopie, 

Universite de Liege, Bat. B15, B - 4000 Liege, Belgium 
4 Instituto de Investigaciones Fisicas de Mar del Plata (1FIMAR), 

CONICET-UNMdP, Funes 3350, B7602AYL Mar del Plata, Argentina. 

5 Consejo Nacional de Investigaciones Cientificas y Tecnoldgicas (CONICET), Argentina 
6 LPTMS, CNRS and Universite Paris-Sud, UMR 8626, Bat. 100, 91405 Orsay, France 

(Dated: June 18, 2015) 

We present a comprehensive study of the destruction of quantum multifractality in the presence of perturba¬ 
tions. We study diverse representative models displaying multifractality, including a pseudointegrable system, 
the Anderson model and a random matrix model. We apply several types of natural perturbations which can 
be relevant for experimental implementations. We construct an analytical theory for certain cases, and per¬ 
form extensive large-scale numerical simulations in other cases. The data are analyzed through refined methods 
including double scaling analysis. Our results confirm the recent conjecture that multifractality breaks down 
following two scenarios. In the first one, multifractality is preserved unchanged below a certain characteristic 
length which decreases with perturbation strength. In the second one, multifractality is affected at all scales 
and disappears uniformly for a strong enough perturbation. Our refined analysis shows that subtle variants of 
these scenarios can be present in certain cases. This study could guide experimental implementations in order 
to observe quantum multifractality in real systems. 

PACS numbers: 05.45.Df, 05.45.Mt, 71.30.+h, 05.40.-a 


I. INTRODUCTION 

Many physical systems display patterns that repeat them¬ 
selves faithfully at every scale. When such systems are char¬ 
acterized by a single non integer dimension, they are called 
fractals (see e.g. yj). More generically, multifractality cor¬ 
responds to the case when different fractal dimensions are re¬ 
quired to describe the system. Multifractality characterizes 
many complex classical phenomena: stock option analysis 
H], turbulence @1, cloud imaging 0]. In quantum physics a 
seminal example where multifractality occurs is the Anderson 
model for the transport of an electron in a disordered crys¬ 
tal 01 . In the metallic phase the electron wave functions are 
spread uniformly inside the sample whereas in the insulator 
phase they are strongly localized. Exactly at the threshold of 
the transition wave functions show highly nontrivial fluctua¬ 
tions leading to anomalous transport. These fluctuations can 
be precisely described by a multifractal analysis, see e.g. 0 
and references therein. Such types of multifractal wave func¬ 
tions can also be found in dynamical systems whose classical 
limit is neither integrable nor fully chaotic, which are dubbed 
pseudointegrable systems ffl.‘ Quantum multifractality in 
various related systems has been intensively studied on the 
theoretical side, from a condensed matter perspective 00- 
0 for both one-body and many-body models and from a 
semiclassical point of view Il20l - [33ll . However experimental 
characterization of multifractality has been much more chal¬ 
lenging, despite some indirect recent attempts in disordered 
conductors 134] and cold atoms 135140 . It is worth mention¬ 
ing that a recent acoustics experiment simulating the Ander¬ 
son model has allowed such a measurement & 

As multifractality has been difficult to observe experimen¬ 


tally, it is crucial to assess how it is affected by perturbations. 
This analysis is also important from a fundamental viewpoint, 
since disturbances of the system may affect the wave function 
at different scales. Considering that multifractality is a multi¬ 
scale phenomenon, this could lead to a wealth of possible be¬ 
haviors. The main goal of the present paper is thus to analyze 
how quantum systems with multifractal properties behave un¬ 
der the effect of an external perturbation. We have considered 
three paradigmatic one-body models, one being the Power law 
Random Banded Matrix model (PRBM), the second one the 
Anderson model, and the third one being representative of 
pseudointegrable systems. In these systems, we have inves¬ 
tigated several natural perturbations in order to specify the ro¬ 
bustness of quantum multifractality. At the same time these 
natural perturbations could account for real experimental situ¬ 
ations. We have recently conjectured that quantum multifrac¬ 
tality can be in general destroyed by a perturbation following 
two scenarios 114111 . In scenario I, there exists a characteristic 
length below which multifractality is unchanged; the pertur¬ 
bation acts only by changing the characteristic length. In sce¬ 
nario II, multifractality is affected at all scales and vanishes 
uniformly when the perturbation increases. In the present pa¬ 
per, we confirm these two broad scenarios by new detailed 
analytical and numerical results. We also introduce a dou¬ 
ble scaling analysis to describe a variant of the second sce¬ 
nario where a modified multifractality is observed only below 
a characteristic scale. 

In Sect.[U]the models we have studied are more precisely in¬ 
troduced and the numerical methods used to obtain our results 
are described. In Sect. UD] we consider a first type of perturba¬ 
tion natural for pseudointegrable models, namely the smooth¬ 
ing of singularities in the potential. In Sect. [IV] we consider 
a change of parameters which moves the system away from 
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criticality. In the case of a specific pseudointegrable system, 
we are able to predict the change of multifractality through an 
analytical theory that we expose in detail. In Sect. [V] we study 
the perturbation corresponding to a change of basis. Eventu¬ 
ally we draw some conclusions in Sect. [VT] 


II. MODELS AND METHODS 
A. Models 

Many theoretical investigations on multifractals were first 
carried out on the example of the PRBM model j 14j] (see also 
Ed). This model is defined (in the real periodic case) as 
the ensemble of symmetric N x N matrices with random real 
coefficients, with zero mean value, and a variance given, for 
1 < i. j < N, by 


{H%) = 1 , = 

(1) 

The parameter b (effective band width) allows to tune the mul¬ 
tifractality of the model from a regime of strong multifractal¬ 
ity (6 < 1 ) to weak multifractality, where states are close to 
extended (b 1). We will use this model as a benchmark 

at specific places, especially since some analytical results are 
available lad]- However, it is not related directly to physical 
models, and in most of the following we will concentrate on 
two other models of more immediate physical relevance. 

The second model we consider originates from semiclassi- 
cal physics | 2 Cj] (see also f' 2 lll) and describes the discrete time 
dynamics of one quantum particle kicked in one dimension 
with a classical limit between integrability and chaos (pseu- 
dointegrability). This model, called the intermediate map, is 
defined as the quantization of an interval-exchange map on the 
torus. The classical map is defined by 

Pn+1 = Pn+ 7 mod 1 , 

x n+ i = x n + 2p n+ i mod 1 . ( 2 ) 

It is generated by the following Hamiltonian, defined on the 
phase space as 

H(p,x) =p 2 + V{x)^26(t-n), (3) 

n 

with V{x) = — 7 {a;}, where {x} means the fractional part of 

x. 

For integrable systems motion in phase space is restricted to 
tori (surfaces of genus one), while for pseudo-integrable sys¬ 
tems motion takes place on surfaces of higher genus. For the 
classical intermediate map with rational 7 = a/b, motion with 
initial momentum po is restricted to the b one-dimensional tori 
(circles) p = po + kj with 0 < k < b — 1 , thus describing a 
surface of genus b. For irrational 7 the motion is ergodic as in 
chaotic systems, although no strong chaos is present. 

The corresponding quantum map is a unitary operator U on 
an TV— dimensional Hilbert space. For the intermediate map it 


1 + 


N . (ir(i — j) 
7T b Sm ( N 


is given in the momentum basis by the N x N unitary matrix 


Uu 


^ — 2Trik 2 /N ^ g2 in'yN 


0 < k,l < N- 1. 


(4) 

The dimension of the Hilbert space is related to the effective 
Planck constant h e g = 2i t/N. We also consider a random 
version of the model, where e~ 2mk ! N is replaced by 
[ l29ll . with 6k independent random variables uniformly dis¬ 
tributed in [0; 27 t]. This model allows to get better statistics, 
and gives similar results as the non-random model, with some 
specificities that we will present. 

The spectral statistics of the quantum map Q depend on the 
value of the parameter 7 . For irrational 7 , the spectral statis¬ 
tics follows the prediction for the Circular Unitary Ensemble 
(CUE) of random matrices characteristic of chaotic systems. 
For rational 7 = a/b, the spectral statistics depend on the 
arithmetical properties of b and are intermediate between the 
Poisson statistics of integrable systems and the random ma¬ 
trix result of chaotic systems 11291130(1 . In the case where 7 
is a rational number 7 = a/b, the eigenvectors of the oper¬ 
ator (|4]i in the momentum basis show multifractal properties 
ll33ll . The multifractality strength depends on b, from strong 
multifractality (small b) to weak multifractality (large b). 

The intermediate map corresponds to the quantization of a 
dynamical system. It is also known that multifractality can 
appear in the critical regime of disordered solid-state systems. 
To discuss this class of systems, we will consider the famous 
model proposed by Anderson in Jit]. The d-dimensional An¬ 
derson model is defined in the basis of lattice sites as 


H = Yl €l l*)(*l + 5Z l*)^l ’ (5) 

* (».j) 


where the random on-site energies e, are uniformly distributed 
in [— W/2, W/2] and (i, j) denotes nearest neighbors. Eigen¬ 
states of this model Q are always exponentially localized in 
dimension one and two. The situation is different in three 
dimensions. Indeed for d = 3 all eigenvectors are local¬ 
ized for large values of the disorder strength W, but the sys¬ 
tem performs a localization-delocalization transition at a value 
W c ~ 16.53 d. For W < W c , eigenstates in the vicinity of 
E = 0 are extended. At the transition point W = W c , states 
display multifractal properties E]. 


B. Multifractal dimensions 

There are several ways of defining multifractal dimensions 
for quantum states, which in most cases yield similar results 
( 3l. In the present paper we will mainly use the box counting 
method. A system of linear size L is decomposed into boxes 
of size t, and a coarse-grained measure of each box k for a 
wave vector \i/>) is defined as /ifc = X^iefc IV’il 2 - We define the 
moments P q (£) of order q as 

= ( 6 ) 
k 
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FIG. 1: (Color online) Multifractal dimensions D q (top left) and sin¬ 
gularity spectrum f(a) (right) for the intermediate map for 7 = 1/3 
and N = 2 12 ; blue solid line: annealed exponent; red dashed line: 
typical. Bottom left: correlation function R 2 {r)\ dashed gray line 
corresponds to the exponent D 2 — 1 obtained from top left. 


In the limit of vanishing ratio t/L between the box size and 
the system size, the presence of multifractality is character¬ 
ized by the following behavior 

^W~(£)\ e/L^ 0 , (7) 

with a nontrivial exponent r q . The main quantity which will 
be used throughout this paper is the multifractal dimension 

D q = Tq/(q~l)- 

Another way of characterizing multifractality is to use the 
scaling of the moments as a function of the system size sin. 
In the systems we study, this method has been shown to be 
equivalent to the box-counting method [33]. Here it will allow 
an analytical approach to be developed, which will be used in 
Section IV B. Nevertheless, this method can be delicate to 
use in certain cases, especially when looking at wave packets 
642]]. Another drawback is that in the systems we consider 
this approach makes it difficult to distinguish the physics at 
different scales, which is crucial in our study. 

Another signature of multifractality is the behavior of cor¬ 
relation functions such as the 2 —point correlation function 

i? 2 (r) = iV 2 (|^| 2 |^+r| 2 ), (8) 

where N = L d is the Hilbert space dimension for a system of 
dimension d, and the average is taken over different eigenvec¬ 
tors, disorder (when present) and all indices i. It is related to 
the multifractal exponent D 2 dlH via 

f? 2 (r) ~ r^ 2-1 , 1-0. (9) 


exponents coincide over a relatively large range of 7 -values in 
the vicinity of q = 0 |33]. As an illustration, Fig.|T|shows an 
example of multifractal dimensions and singularity spectrum 
for the intermediate map. In what follows we will mainly con¬ 
centrate on the set of annealed exponents. We also display in 
Fig. |T|the correlation function f ? 2 for the intermediate map, 
together with the slope corresponding to Z) 2 , illustrating rela¬ 
tion (|9}. 


C. Local multifractal exponents 

Multifractality is mathematically defined as a scale invari¬ 
ance which takes place at all scales. In a real setting however, 
multifractality can be valid only on a certain limited range of 
scales, e.g. between a lower microscopic length and an upper 
macroscopic length. As we will show, this is particularly rel¬ 
evant for perturbed systems. In order to investigate the ways 
in which multifractality is destroyed when a system is per¬ 
turbed, one can introduce Ilia i4ltl a local multifractal expo¬ 
nent, which characterizes multifractality at a given scale. It is 
defined as 


D q {t) 


1 din P q {t) 
q — 1 d In £ 


( 10 ) 


In practice as the scales £ are discrete numbers we compute 
the local multifractal exponent at scale £ as the slope between 
scale £ and the scale immediately above. 

The local multifractal exponents typically show a plateau 
which corresponds to the global multifractal exponent defined 
by (0: as an example. Fig. [2] displays the local multifractal 
exponents for the PRBM model, together with analytical pre¬ 
dictions from [60, and Fig.[3]displays this quantity for the in- 



£ 


Alternatively, one may express multifractal properties via 
the singularity spectrum /(a), which is the Legendre trans¬ 
form of T q . Moreover, for disordered systems, one distin¬ 
guishes between the annealed exponents which describe the 
scaling of the average moments, and the typical exponents 
which characterize the average of the logarithm of the mo¬ 
ments. For all the systems considered here, the two sets of 


FIG. 2: (Color online) Local multifractal dimension D 2 (£) as a func¬ 
tion of the box size £ for the PRBM model CD for two values of 6. 
Top: weak multifractality regime 6 = 2. Bottom: strong multi¬ 
fractality regime 6 = 0.01. Black circles: N = 2 10 ; red squares: 
N = 2 11 ; green diamonds: N = 2 12 ; blue triangles: N = 2 13 . 
Dotted dashed orange line: analytical prediction D 2 = 26 for 6 <C 1 
and D 2 = 1 — 75 for 6 » 1 0. 
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FIG. 3: (Color online) Local multifractal dimensions for the inter¬ 
mediate map for 7 = 1/5 and increasing system sizes N = L. The 
full lines correspond to q = 2 while the dashed lines correspond to 
<7=1. Black circles: N = 2 9 ; red squares: N = 2 10 ; green dia¬ 
monds: N = 2 11 ; blue up triangles: N = 2 12 ; orange left triangles: 
N = 2 13 ; brown down triangles: N = 2 14 . Left: raw data. Right: 
all the box sizes are rescaled by S = N/b. 

termediate map. Both figures show that D q {() indeed presents 
a plateau for a significant range of £ values. However, devia¬ 
tions can occur for the smallest values of £ as in Fig. [2] which 
are due to the fact that coarse graining is necessary to obtain 
converged results. 

Deviations can also occur at large scales, as is the case for 
the intermediate map in Fig. [3] The plateau at small £ coin¬ 
cides with the value D q from Fig. [7] but at large scales D q {£) 
saturates to 1. This is a specificity of the model, which for 
7 = a/b exhibits a characteristic length 3 = N/b, arising 
from the existence of the underlying classical structure de¬ 
scribed in section III Al The characteristic length can also be 
seen on the data shown on Fig.Q]for the correlation function, 
where b peaks of typical width N/b are clearly visible. Be¬ 
low the characteristic length 3, the value of D q (£) shows a 
plateau indicating asymptotic multifractal behavior. Figure 0 
shows that after rescaling the box size £ by this characteristic 
length, all D q {£) collapse to a single curve following the law 
D q (£) = -F g (£/3), where T q (u) is a function independent of 
N. 


D. Natural perturbations and scenarios 

In the following Sections, several types of natural perturba¬ 
tions will be applied to these different models. For the inter¬ 
mediate map © three types of perturbations naturally arise: i) 
the singular potential in 0 can be smoothed; ii) the parameter 
7 which controls multifractality can be varied away from its 
critical values; and iii) the measurement basis can be changed. 
Clearly, all these perturbations can arise in a real experimental 
setting, such as cold atom lf35l - [39ll or photonic lattice iQ im¬ 
plementations. In the case of the Anderson model (0, the nat¬ 


ural perturbations correspond to a change of disorder strength 
away from criticality and a change of measurement basis. 

We will develop a scaling analysis of numerical data com¬ 
bined with analytical approaches in order to show that the dif¬ 
ferent paths to multifractality breakdown always follow one of 
the two scenarios presented in [41 ] and outlined in the intro¬ 
duction: scenario I corresponds to the existence of a character¬ 
istic length below which multifractality is unchanged; above 
this characteristic length which decreases with increasing per¬ 
turbation strength, multifractality is destroyed. Scenario II 
corresponds to a multifractality which is affected at all scales 
and vanishes uniformly when the perturbation increases. We 
will see that there can be interesting variations depending on 
the interplay between characteristic lengths of the model and 
the perturbations. 


III. SMOOTHING THE SINGULAR POTENTIAL 

In this section several types of smoothing of the intermedi¬ 
ate map are described, which aim to account more realistically 
for experimental constraints. Indeed, in pseudo-integrable 
systems, generally singularities are present and are one of the 
reasons for which the classical dynamics is neither integrable 
nor chaotic. Experimentally, discontinuities such as in the po¬ 
tential in 0 have to be smoothed out. In this section we will 
consider several types of smooth potentials approximating the 
exact one in different ways. These perturbations have been 
thought to be relevant for an experimental implementation of 
the intermediate map. One could envision photonic crystal 
implementations |4J] where time is taken as a spatial dimen¬ 
sion and the potential is etched on a substrate whose refrac¬ 
tive index is varied. In this context, the potential singularity 
will be smoothed over a certain distance which depends on 
the etching technique. Another possible implementation cor¬ 
responds to cold atom experiments where atoms are subjected 
to potentials constructed from laser light standing waves |35- 
13911 . In this context, the smoothing of the potential will take 
place through the presence of only a fraction of the Fourier 
components needed to build the exact potential in 0. Three 
possible ways of smoothing the potential V (x) are considered 
below, adapted to these two experimental possibilities. We 
consider both the model with random phases and the deter¬ 
ministic model 0 (see section llTAl) . more realistic for exper¬ 
iments. 


A. Polynomial smoothing 

We first consider a more realistic version of the model for 
photonics experiments [44]. In this context, we chose to ap¬ 
proximate the potential V (x) as 

V(x) = [ ~ 1X ’ 0 < x < 1 - e 

^ ' y a^x^ + CL 2 X 2 + a\X + ao, 1 — e < x < 1 ’ 

( 11 ) 

where the a* are chosen to make the potential and its first 
derivative continuous at x = 1 — e and x = 1. The origi- 
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FIG. 4: (Color online) Intermediate map smoothed by a polynomial 
for (V = 3 9 and 7 = 1/5. Top: Smoothed potential. Dark brown 
full line: exact potential. Dashed purple line: e = 0.05. Dot dashed 
magenta line: e = 0.1. Bottom: Moment of the eigenvectors of 
the random intermediate map for q = 2 as a function of box size. 
Dark brown full line: exact potential; dashed green line: t = 0.005; 
dashed dotted purple line: e = 0.05. 


nal model (0 is recovered when e = 0 so that e can be seen 
as a small perturbative parameter. Typical examples of the 
resulting potential are shown in Fig. [4] top, while typical re¬ 
sults for the moments P q {£) of the random intermediate map 
are shown in Fig. [4] bottom for N = 3 9 . In Fig. [5] the local 
multifractal exponent L >2 (£) is plotted as a function of £ for 
several smoothing widths e. In the left panel the raw data are 
shown: while at very small perturbation strength e one ob¬ 
serves a plateau at the unperturbed D 2 (compare with Fig. [3j. 
at larger values of e this plateau is no longer visible and the 
curves l) 2 {£) increase monotonic ally. Nevertheless, it turns 
out that one can put all these different curves onto a single 
one by rescaling the lengths £ (right panel). This shows that 
the local multifractal exponents obey a scaling relation: 


D q {£) = g q , ( 12 ) 



FIG. 5: (Color online) Scaling analysis of the local multifractal di¬ 
mension D 2 for the perturbed random intermediate map for N = 3 9 
and 7 = 1/5. The perturbation is the smoothing of the singular 
potential over a length e. Left: raw data. Right: data after rescal¬ 
ing. Inset: variation of the scaling length £ as a function of e; solid 
line is the fit corresponding to GU with a = 1.04. Black circles: 
e = 0.002; red squares: e = 0.003; green diamonds: e = 0.005; 
blue up triangles: e = 0.008; yellow left triangles: e = 0.01; brown 
down triangles: e = 0.02; gray right triangles: e = 0.03; purple 
plus: e = 0.05; cyan crosses: e = 0.07; magenta stars: e = 0.1. 

with £ (e) a scaling length which depends only on the pertur¬ 
bation strength and which is well fitted by 

£(e) « e- Q (13) 

with a w 1 (see inset of Fig. 0 and Q q is a scaling function 
independent of e. This scaling behavior is valid for various 
values of the parameters q and N. Indeed, it was shown in 
El to occur when N is of the form 2" while Fig. 0 shows 
that it remains also valid for N of the form 3". Moreover, this 
scaling behavior with exponent a « 1 applies for different 
values of q and 7, and also for polynomial smoothings CD of 
higher order (data not shown). 

Multifractality in this perturbed model can be further inves¬ 
tigated using the 2 —point correlation function Pair) as de¬ 
fined in (0. In Fig.0 is shown for several smoothing 

widths. Rescaling r by the same scaling parameter £(e) oc 1/e 
leads to the collapse of all the curves, see Fig.0right. 

A similar behavior can be observed for the model (0 with 
non-random phases. In Fig.0the local multifractal exponent 
D 2 (£) for the deterministic model is shown to follow the scal¬ 
ing law (IT 2 l >. in complete analogy with Fig.0 

The noticeable difference lies in the exponent a of the scal¬ 
ing length £(e) with respect to the smoothing width e, which 
turns out to be a = 0.67 rather than « 1 for the random phase 
model. As in the case of the random model, this value of 
a ~ 0.67 does not depend on N or 7 . Results for the 2—point 
correlation function are shown in Fig. 0 using the same pa¬ 
rameters as in Fig.0 Again the same scaling behavior is ob¬ 
served, with an exponent a = 0.67 for the rescaling of r. The 
difference of the scaling exponent a between the random and 
non-random model reflect the differences of the correlations 
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FIG. 6 : (Color online) Two-point correlation function R 2 as defined 
in © for the random intermediate map smoothed by a polynomial 
for N = 3 9 and 7 = 1/5. Left: raw data. Right: data after rescaling 
the variable r by £ = 1/e following 113b . The same color code as in 
Fig.[5]i s used. 


FIG. 8 : (Color online) Same as Fig. [ 6 ] for the non-random inter¬ 
mediate map for N = 3 9 and 7 = 1/5. Left: 2—point correla¬ 
tion function as defined in ©. Right: the variable r is rescaled by 
£(e) = e -0 ' 6 ', different from the scaling used in Fig. [ 6 ] The same 
color code as in Fig.[5]is used. 



FIG. 7: (Color online) Same as Fig.[5]for the non-random intermedi¬ 
ate map for (V = 3 9 and 7 = 1/5. Left: raw data of the local mul¬ 
tifractal dimension D 2 for different perturbation strengths. Right: 
data after finite-size scaling. Inset: Variation of the scaling length 
£ as a function of e; solid line is the fit corresponding to O with 
a = 0.67, different from the result in Fig. [5] The same color code as 
in Fig.[5]is used. 

in the phases of the propagator coefficient ©. 

The data discussed in this section show that this kind of 
smoothing leads to the appearance of a characteristic length 
below which multifractality is unchanged, indicating that in 
this case multifractality breakdown occurs following scenario 
I of Subsec tion lll Dl 


trolled with a very high accuracy. One could think that each 
Fourier component of a periodic potential can then be simu¬ 
lated by one laser so that any potential could be reproduced. 
The limitation is that it is practically impossible to use a large 
number of lasers so that only potentials with a small number 
of non zero Fourier components can be modeled. The poten¬ 
tial of the intermediate map is acting on the torus so it can be 
expanded as a discrete Fourier series. In this section we will 
investigate how the eigenvector statistics changes when the 
Fourier expansion of the potential is truncated to Nf terms. 

For Nf = N the linear form is recovered. A plot of the 
potential for different values of Nf is shown in Fig. [9] Con¬ 
trary to the preceding case, the modification of the potential 
is not local anymore. In particular, even for large values of 
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B. Fourier series smoothing 

In a cold atom experiment the potential experienced by the 
atoms can be created with standing waves of laser light. In 
these setups the frequencies and the amplitude can be con- 


FIG. 9: (Color online) Potential of the intermediate map obtained by 
truncation of Fourier series for 7 = 1/5 and N = 3 s = 6561 for 
different Nf. Black full line: Nf = N = 6561; red squares: Nf = 
6559; green diamonds: Nf = 6557; blue up triangles: Nf = 6555; 
orange left triangles: Nf = 6001. Inset is a blow up of the same 
data. 
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FIG. 10: (Color online) Two-point correlation function for the ran¬ 
dom intermediate map with potential approximated by a truncated 
Fourier series for N = 3 8 and 7 = 1/5. Contrary to Fig. [6] we 
do not observe a systematic dependence on the truncation, indicating 
that the approximation is never perturbative. The same color code is 
used as in Fig. [9] 

Nf oscillations remain visible far from the discontinuity, see 
the inset in Fig. [9] In this case our investigation shows that 
even for Nf close to N, when almost all the Fourier compo¬ 
nents are kept, multifractality is completely destroyed. This 
indicates that this kind of perturbation is always large, and 
cannot be made arbitrarily small due to the discreteness of the 
Fourier series. This is illustrated by considering the 2-point 
correlation function in Fig.|T0j contrary to Fig. [6] we do not 
observe a systematic dependence on the perturbation strength, 
even for Nf close to N. Our results therefore show that a 
naive truncation of the Fourier series is not a good approach 
to experimentally observe multifractality in such systems. 


FIG. 11: (Color online) Potential of the intermediate map approx¬ 
imated by a trigonometric smoothing with Nd derivatives fixed at 
x = 0.5, N = 2 12 and 7 = 1/5. From right to left on the left: 
Nd = 46 (red), Nd = 80 (green), Nd = 109 (blue); black straight 
line is the exact potential. 
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C. Trigonometric smoothing 

In view of the results of the preceding subsection, one may 
try to search for better approximation schemes using a mod¬ 
ified Fourier expansion. Indeed, a more efficient way to ap¬ 
proximate the potential for cold atom experiments can be de¬ 
vised by fixing a prescribed number N r j of derivatives at one 
point in order to force the potential to be approximately lin¬ 
ear around that point. Such a potential can be chosen as a 
trigonometric sum: 

K 

V{x) = £ ai sm.{Ttlx) . (14) 

i=o 

We took K = 3 A// and the a/ are fixed by the Nd equations 

y'(0.5) = -7, V"(0.5) =0, ,..., y (JVd) (0.5) = 0 . (15) 

The resulting potential for several values of Nd is shown in 
Fig. HU Compared with the simple truncation of the Fourier 
series, the potential change now occurs only in a limited re¬ 
gion of space, which gets smaller and smaller as Nd increases. 


FIG. 12: (Color online) Scaling analysis of the 2—point correlation 
function for the random intermediate map with potential approxi¬ 
mated by a trigonometric smoothing with Nd derivatives fixed at 
x = 0.5 for N = 2 12 and 7 = 1/5. Black full line: exact model. 
Red circles: Nd = 46. Green squares: Nd = 80. Blue diamonds: 
Nd = 109. Left: raw data; right: data rescaled with dl6l with 
e = 1 /N d . 


We found that a scaling analysis similar to the one in Subsec - 
tion llll Al is possible in this case. As an example, the two-point 
correlation function for different A// is displayed in Fig. [12] 
The scaling relation follows the formula 

R 2 (r) = H , £(e) oc i , (16) 

where e = 1/Nd and TZ is a scaling function independent of 
Nd- This is similar to what is described in Fig. [6] 

This way of expanding the potential as a trigonometric se¬ 
ries is thus more efficient in order to keep the multifractality 
of the system, and the disappearance of scale invariance cor¬ 
responds to scenario I (see Subsection lllDk 
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In all the models considered, multifractality is predicted 
only for certain critical values of a parameter. In this Sec¬ 
tion, we investigate the robustness of multifractal properties 
when this parameter is moved away from criticality. 


A. Change of W in the Anderson model 

In the case of the 3D Anderson model, multifractality ap¬ 
pears at the metal-insulator transition which corresponds to a 
specific disorder strength W c ~ 16.53 d. A natural choice 
of perturbation is therefore to change the disorder strength 
slightly below or above the critical value W c . In this case, it is 
known that the eigenstates are either localized or delocalized 
with a characteristic length £. In the insulating phase £ corre¬ 
sponds to the localization length, while in the metallic phase 
it corresponds to the correlation length. Below this character¬ 
istic length, the wave functions are multifractal with the same 
critical multifractal spectrum, and they form a “multifractal 
insulator” or a “multifractal metal” 1450 . This is a conse¬ 
quence of a one-parameter scaling law that governs the mul¬ 
tifractal spectrum in the vicinity of the transition jjjf]. This 
type of perturbation therefore follows scenario I of quantum 
multifractality breakdown: quantum multifractality survives 
unchanged below a certain characteristic length related to the 
distance to the critical point. 


B. Change of 7 in the intermediate map 

1. Numerical results 

The behavior of the intermediate map is richer, in the sense 
that there is in principle an infinite number of critical values of 
the parameter, i. e. values for which multifractality arises. In¬ 
deed, in the intermediate map i(4]>. multifractality is predicted 
to appear for rational values of the parameter 7 (see Section 
HB . Multifractality manifests itself all the more strongly when 
the denominator of 7 is small, that is, for 7 = 1/2,1/3 or 
1/5 for instance 1131 11 . In close vicinity of these rationals, 
one should observe delocalized eigenstates. However, since 
this difference in behavior only arises in the limit of infinite 
size, we can expect at finite size N a persistence of multi¬ 
fractal properties if one varies the parameter in some vicin¬ 
ity of these low-denominator rationals. This is illustrated in 
Fig. Qj] where D q is computed for a fixed vector size N as 
a function of the parameter 7 in the vicinity of two rationals, 
1/2 and 1/5, up to a distance of the order 1/N from these 
rationals. Clearly the curve D q ( 7 ) is not singular at all ra¬ 
tionals, but rather smoothed out. An advantage of this model 
is that it is amenable to analytical treatment via perturbation 
theory, which allows us to get a clear picture of how parame¬ 
ter changes may affect multifractality. This approach will be 
carried out in the next subsection. 
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FIG. 13: (Color online) Variation of the multifractal dimensions 
D± 1 ( 7 ) for the intermediate map (01 with random phases and N = 
2 12 in the vicinity of the rational values 7 = 1/2 (top) and 7 = 1/5 
(bottom). The red solid line below 1 corresponds to D 1 ( 7 ), the blue 
one above 1 to D_ 1 ( 7 ). Black dash-dotted parabolas correspond to 
the theoretical expression Eq. 0 D> with k = 1 and 2 for 7 = 1/2 
and k = 0 and 1 for 7 = 1/5. Inset: zoom on the right part of the 
top plot, corresponding to k = 2 . 


A crucial property of this type of perturbation, as our nu¬ 
merical results show, is that multifractal properties for differ¬ 
ent values of 7 do not depend on any characteristic length. 
This can be shown by investigating the 2—point correlation 
function © for this model. Results displayed in Fig. |T4] ftop) 
show that f ?2 behaves as a power law as in Eq. © over a broad 
range of scales and yields a well-defined multifractal dimen¬ 
sion £>2 which increases toward the ergodic value £>2 = 1 
when the parameter is tuned away from criticality. The fact 
that multifractal dimensions change smoothly and uniformly 
at all scales is a footprint of scenario II. To confirm this ef¬ 
fect we have also computed higher order correlation functions, 
which are known to involve other multifractal dimensions @|. 
An example is shown in Fig.fl4l(bottom). showing that indeed 
other multifractal dimensions follow scenario II. 
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FIG. 14: (Color online) Correlation functions for the random inter¬ 
mediate map of size N = 2 12 in the vicinity of 7 = 1/5. Top: 
2 -point function §}. bottom: higher order correlation function. The 
black dashed straight line corresponds to r 2D3_D2_1 , as expected 
from ||4j]. The curves correspond to 7 = 1/5 + e/(5N) with (from 
top to bottom) e = 0 (black solid line), 0.127 (blue dashed line), 
0.255 (green dotted line), 0.382 (orange dash-dotted line), 0.509 
(dark-red dash-dotted line), 0.637 (light-blue double-dashed line), 
0.764 (purple dash-double dotted line), 0.955 (red dashed line). 


2. Outline of the perturbation theory 


As mentioned above, it is possible to build analytically a 
perturbation theory to describe the vicinity of rational values 
of 7. This is done by using a related mathematical_model, 
the Ruijsenaars-Schneider model [[46]. As shown in [32], this 
model also displays multifractal eigenvectors. It was already 
used to describe a quantum version of the intermediate map 
with unbounded phase space [42] - 

In this section we develop a perturbation theory for the in¬ 
termediate map when the parameter 7 depends on the matrix 
size. Namely, we consider the intermediate map at parameter 
value of the form 7 = 1/6 + a/N, that is, when the parameter 
of the map gets close to a rational at a speed which depends 
on N. The Ruijsenaars-Schneider (RS) ensemble is defined 


as the ensemble of N x N unitary matrices of the form 

_ gifm 1 _ e 2mg 

Umn — 1 _ e 27 ri (m-n+g)/N ’ 

with independent random phases uniformly distributed 
between 0 and 27 t, and g some fixed parameter 113211 . 

The perturbation expansion for multifractal dimensions of 
the Ruijsenaars model was obtained in [[32;!] in the weak mul- 
tifractality limit where g is close to a nonzero integer. The in¬ 
termediate map at parameter value 7 = 1/6 + a/N coincides 
with the Ruijsenaars map with parameter g = N 7 = N/h+a. 
Let r be the remainder of N modulo 6. The weak multifrac- 
tality limit for ( fI7l ) is obtained when N/b + a is close to a 
nonzero integer, that is, when a = k — r/b + e, with k an 
integer and e a small real number. For such a value of a, the 
intermediate map corresponds to the map ( IT71 ) with parameter 
g = k + e, where R = (N — r)/b + n is an integer. Thus we 
expect multifractal dimensions for the intermediate map to be 
given by a perturbation expansion in e similar as in [32], but 
in the vicinity of an integer k which depends on N. For sim¬ 
plicity, we will consider the case where gcd(K, N) = 1. We 
can then define k - 1 as the inverse of k modulo N. 

Multifractal dimensions can be obtained from the asymp¬ 
totic behavior of 

^2 \^ri{0t)\ 2q ~ ( 18 ) 


at large N. Here T„ (o:) is the nth component of the nth 
eigenvector of the system ( ITTb of size N. In the unperturbed 
case where states are extended, the multifractal dimensions 
are I)',/ 1 = 1 for all q. The perturbative approach allows to 
express eigenvectors of ( fT7] i at first order in e, and thus the 
moments of the wavefunction. At e > 0, fractal dimensions 
are given by D q = 1 — \d q , with d q some small number (the 
factor 7j is put here for convenience). From (fl8l) one then ob¬ 
tains 


^|'L„(a)| 29 ~ N -0-Qd g /2)( q -i) 
n 

The first-order correction to the multifractal dimension is thus 
given by the logarithmic behavior of the perturbative correc¬ 
tion of the moments. We first find a closed expression for 
the first-order correction of the wavefunction moments aver¬ 
aged over the whole spectrum and over disorder configura¬ 
tions (Eqs. ( |28| > and 07]). and then extract the dominant log¬ 
arithmic contribution in the limit N —> 00 (Eq. (l42l >). which 
gives us the correction d q sought for. 

Eventually, the calculations detailed in the next paragraph pro¬ 
vide an analytical confirmation that a change of 7 in the inter¬ 
mediate map leads to a multifractality breakdown following 
scenario II. 
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3. Perturbation expansion 

Let us consider a perturbation expansion of (ITTb around 
k, setting g = k + e. Let M mn = U mn e~ lne( ' 1 ~ 1 ^ N ' l \ this 
rescales U mn by a trivial factor, and 


where terms linear in q sum up to zero because of the normal¬ 
ization of Identifying ( l28l ) and ( IT9l ). we get 

jj^iQnia) + Q* n (a)) 2 ~ d q ]nN, (29) 


Mmn — &m-n-\-K, ~ 


■ 

e m sm7re 


4 “ (1 fim—n+k) 


N sin(7re/lV) 

_ g2Trie^ e -i7r€(l-l/tV) 
g27 ri(m—n-\-k-\-e)/N 


( 20 ) 


is then such that both terms have a definite limit when e —> 0 . 
First-order expansion of M mn reads 


JT “ 1 " 1 v 

mn — & Om-n-\-k 77”^ 


2 z 7 T€ i$> m 1 $m-n-\-k 

n+k N ~ 6 l — e 2 iri(m-n+k)/N ' 

( 21 ) 

We denote eigenfunctions and eigenvalues of M mn respec¬ 
tively by v F n (a) and A a . Unperturbed eigenstates, that is, 
eigenvectors of e l4>m are given by 


^ 0) («) = 


-S„(a) = —to + ^ $kj 

1=0 


with eigenvalues 


A (0) = gii+^a 


( 22 ) 


(23) 


1 \-^N — l 
N 

expansion of eigenvectors gives 


where $ = jt Ej=o Standard first-order perturbation 


*»(«) = *£>(«) + ]T C a ^\p), (24) 


so that the correction to the unperturbed multifractal dimen¬ 
sion is indeed given by the logarithmic asymptotic behavior 
of the sum in 1 


4. Model with random phases 
From the definition d27l > of Q n we have 

“^7 ^ ' Qnipt) = ^ ' Ca,ot-\-filCa,a—l3 ; 


(30) 


and 


N 


J2 Qn(a)Q*Ja) = C a , a+ pC* a>a+fi . (31) 


We are interested in quantities averaged over all eigenvectors 
and random phases: we thus have to perform a sum over a 
and an integral over the <1> 7 . From Eqs. ( 12 1 1 ( 1251 ). the explicit 
expression of C a> p+ a is 


C a ,/3+a = 2lV £r S eX P *( r ~ s ~ 1 )(^F + $) 


x exp 


-^f( s + \)P - i Ej=o + * Ej=o 


(32) 


(we have changed the summation from m, n to r, s with r = 
k~ 1 n and s = k~ 1 m), and 


with 


'-'a/3 — 




\(0) _ x (0) 


(25) 


and Mmn the order-e (off-diagonal) term in (f2Tb . Replacing 
by its explicit value d22l > in Eq. (l2-fl > we get 

'I'n( a )| 2 = 177 (1 + Qn{°) + Qn ( a ) + Qn{&)Qn{ a )) 

(26) 


N 


with 


Q n {a) = Y, e2 ™ MNc ^+° 


(27) 


Taking (1261 ) to the power q and summing over n we get, up to 
order e 2 , 

E r■ E («»(“) + «(«)) 2 . 


ne e -inx/N 

t x = — r —- 7 — if x ^ 0, 0 otherwise. (33) 

N santx/N 

Similarly C a}0l -p can be expressed by a sum over indices r' 
and s'. The quantity C at p-\- a C at0l -p depends on a through a 
factor exp[j(r — s — 1 + r' — s' — l)^r]- Upon averaging 
over a, one thus gets a coefficient 6 r - s + r '- s '-2 which kills 
the terms <h. The averaging of C 0 n:,+HC n m-a over random 
angles then contains a coefficient 


exp 


r' — 1 


1=0 1 =0 1=0 1=0 


(34) 

Since r — s + r' — s' — 2 = 0, and to = 0 (so that contributions 
with r = s + 1 or r' = s' + 1 vanish), the average (l34l) can 
only be nonzero when s = r' — 1 and s' = r — 1. This yields 


C, 


a , a+/3 C a , ck—/3 


_ _ l^(s —r~+l)K I 2 | 2 L( s _ r+1 )g 

41V 2 ^ ^ 


( 28 ) 


( 35 ) 
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Changing variables s — r + 1 = x and summing over /?, we 
get from (l30l > and (l35l > 


/ 1 _ 7T 2 e 2 1 e 2 n x0 

\ N ^ ^ n[a) /a,* “ 41V 3 sin 2 2*r gin 2 ‘ 

n x,fi N N 

(36) 

In a similar way, averaging C' ajQ , +j gC'* iQ , +/3 over a yields a 
coefficient <5 r _ s _ r / +s /, and the average over then yields 
the condition r = r' and s = s'. An expression equivalent 
to (l36l > can be found, which, summed together with Eq. ( 1361 ). 
yields 


/ 1 


( 1 y yi(Qn(oQ + Qn( a ))‘ 


E 


x(N — x) 


t> IV 3 ^ sin 2 (7rra/)V)’ 
(37) 

where we have performed the sum over j3 by using the identity 


y- 1 sin 2 (n(3x/N) 
sin 2 (n/3/N) 


x(N — x). 


(38) 


From (l29l) we know that the lowest-order contribution to the 
multifractal exponent is given by the logarithmic behavior of 
the sum (l37l > for large N. Recall that k = (N — r)/b + k with 
k fixed. The sum in (l37l > can be split into 6 subsums, as 


2 2 1 IW/ b J t , | w A t , \ 

Y' Y (bx + c){N - bx - c) 

-2 (c i (Kb-r)(bx-\-c)\ 

c=o x=0 sm tt + -- Yb - -) 

(we obviously omit the case c = x = 0 in the above sum). The 
logarithmic contribution ( l29l ) originates from regions where 
the divergence of the sin 2 in the denominator is compensated 
by a linearly vanishing numerator. This corresponds to the 
two regions bx + c ~ 0 and bx + c ~ N (in all other cases, 
either the summand is ~ \/x 2 and converges, or it gives 
non-logarithmic divergences which should be compensated by 
higher-order terms in the perturbation expansion if we assume 
that the behavior (l29l > holds). The first region comes from the 
sum with c = 0 in ( l39l >. the second one from the sum with 
c = r. 

The sum for c = 0 in ( l39l > can be rewritten as the sum of 
two terms, namely 



bN 3 


(nb — r) 2 ir 2 x 


(40) 

which is a Riemann sum converging to an integral with a finite 
value, and 


be 2 

(nb — r) 2 


LJV/bj-i 

E 


X — l 


1 

x 


be 2 

(nb — r) 2 


In N 


(41) 


which is responsible for the logarithmic divergence. After a 
change of variables x (N — r)/b — x, the sum for c = r 
in ( [39b can be shown to yield exactly the same contribution as 
c = 0. Summing both contributions, one gets 


(^E^» + 3») 2 } 

n 


a,<3? 


2be 2 

(nb — r) 2 


In TV, (42) 



FIG. 15: (Color online) Multifractal dimensions D 1(7) for the inter¬ 
mediate map with random phases (blue solid line) and deterministic 
phases (red dashed line) in the vicinity of 7 = 1/2 for N = 2 12 . 
Black dash-dotted parabolas correspond to the theoretical expression 
Eq. (S3} for k = 1, 2 while black dotted parabolas correspond to 
Eq.03 with an additional prefactor 2 in front of q. 


which, by identification with ( l29l ). gives d q and thus D q . Since 

7 = 1 /b + (k — r/b + e)/N, one finally gets 


D q 



(l-N 


76 — 1 
Kb — r 


2 


(43) 


The first maxima of D q around 7 = 1/6 correspond to k 
0, ±1, ±2 ... The dependence in N in Eq. ( l43l > indicates that 
this fractal dimension D q gives the behavior of wavefunctions 
of size N when the parameter 7 is considered at a scale 1/N 
around the rational 1/6. 

The theory can be compared with the numerical results; 
Fig. [13] shows that indeed it describes correctly the vicinity 
of rational points for finite N for positive or negative q. 


5. The case of deterministic phases 

It is instructive to also analyze the fractal dimensions in the 
case where random phases are replaced by the phases <17 = 
2"k 2 /N of the deterministic intermediate model. Numeri¬ 
cally, this model yields results which are close but slightly 
different from those of the model with random phases. While 
the difference is tiny for most values of the denominator 6 of 
the parameter, the most significant discrepancy occurs in the 
case 7 = 1/2. For instance, as illustrated in Fig.[l5j when 7 is 
equal to 1/2 and N is a power of 2, the value of 1 — l) \ for de¬ 
terministic phases is twice the value for random phases given 
by Eq. ( l43l ). The method used above can in fact be adapted to 
explain these discrepancies. In this subsection we will show 
this for 7 = 1/2 and N = 2 n . Other cases can be treated by a 
similar approach. 

To obtain the result for the deterministic model, the main 
change in the perturbative calculation of D q comes from 
Eq. (l34l >. where the average over random phases is replaced 
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FIG. 16: (Color online) Left: Function /n(x) defined in d46b (see 
text) as a function of x for N = 2 11 (black), and function x(N — x) 
(dashed/red). Right: Same plot with a different ordering of the ab¬ 
scissas. Points with abscissa x = 2 n ~ q (2y+l) withg = 1,2,... ,n 
and y = 0,1,..., 2 q ~ 2 are ordered by increasing q and at fixed q 
by increasing y. The last family (from 512 to 1024) corresponds to 
q = n, that is, odd x. 


by a constant term £ rs £ r ' s ', with 

- 2 / 1 l 


£ rs = exp 


2l7TKT 

N 


V k=0 k —0 


(44) 


Performing the same steps as previously, one can show that 


/ 1 


T7 ^2(Qn(a) + Qn(a))‘ 


\N 


7t 2 e 2 /jvQr) 

N 3 ^ sin 2 (nkx/Ny 


(45) 

where the numerator x(N — x ) of Eq. (l37l) is now replaced by 
a function 


labels x by setting x = 2 n ~ q (2y + 1) with q = 1,2,..., n 
and y = 0,1,..., 2 9-2 , we obtain the right panel in Fig. [16] 
We see that the linear behavior of /n{x) does not occur only 
for x ~ 0, but also for all families x = 2 n ~ q (2y + 1) when 
y ~ 0. Namely, for all q and y small we observe numerically 
that 


f N (x) ~ x(N -x) = 2^ n ~ q \2y + 1)(2« - 2y - 1). (48) 


For h = N/2 + 1 and even x, i.e. families x = 2 n q [2y + 1) 
with 1 < q < n — 1, we have 


. 2 TTKX . 2 

sin —— = sm 
N 




= sin 2 


n(2y+l) 

29 


(49) 


so that (]48l > will partly compensate ( l49l ). and thus for small y 
each q-family with q ^ n will give a logarithmic contribution 
to (l45l >. As before, odd x do not contribute, since for x = 
2y + 1 (the family q = n) we get 


. ,7 TKX o TTX n7r(2l/+l) 

sm —— = cos — = cos ---. 

N N 2 n 


(50) 


so that vanishing of <[50t does not correspond to a vanishing 
of d48l >. Using d48l) and d49t . the contribution to d45l of a q- 
family x = 2 n ~ q [2y + 1) with 1 < q < n — 1 for small y 
reads 


^ 2 e 2 2 2 (n- q ) (2y + 1)(29 - 2y - 1) 6 2 ^ 1 

N 3 2^/ n(2y+l) 0 „ 4 


N 3 -2 7T(2j/+l) 

y &m 29 

at 

V 

2y + l 



(51) 

Then for large q 



2 q -2 2«- 1 +l 1 27— 2 

y 1 = y y 1 

^ 0 2 t/ + l ^ y ^2 y 


(52) 


and the sum over all contributions ( |5TT > gives 


fN(x) = ^ cr(s' - s)<5 r _ s _i_ x (46) 

rr' ss' 

^ [^r,s( £,r' ,s' &r'—s' — l-\-x H - £,r',s' &r'—s' — 1 —x)] 5 

with a an N- periodic function defined for 0 < v < N — lby 

N 2 - 1 


= Y, 7: 




y sin (tt/3/N) 


— 2v(N — v). (47) 


As before, the logarithmic behavior of (l45l) is expected to 
come from places where sin 2 (7 tkx/N) becomes close to zero 
while the numerator approaches zero linearly. Previously, 
when the numerator was given by x(N — x), this only oc¬ 
curred for x cs 0 and x ~ N. Here the function /at (a;) 
still has a linear behavior in the vicinity of x ~ 0 but it is 
much more oscillating for larger x. For illustration, we plot 
an example of Jn{x) in Fig. [I6](left panel) for N = 2 n ; we 
concentrate on indices x £ [0, iV/2] since by symmetry the 
other half yields the same contribution. The only linear con¬ 
tribution would seem to come from x ~ 0. However, let us 
consider the case N = 2 n and 7 = 1/2. If we rearrange the 


2 71—1 2 

2 - V2«ln2«~ i-lnJV. (53) 

2 AT 4—( 2 

q= 1 

So far we considered only the region x < N/2. The region 
x £ [N/2,N] contributes in the same way, so that the total 
contribution from the oscillating part of /n{x ) is twice the 
result (f53l) . that is, e 2 In N. To this contribution one must add 
the contribution from the linear part of /n{x), corresponding 
to x ~ 0 (and x ~ N), see Fig. [T6] left; the calculation is the 
same as for random phases and thus yields the contribution 
given by Eq. (l42l >. For b = 2 and k = 1 this term is equal to 
e 2 In N. The total of all contributions for deterministic phases 
is thus 2f 2 In N, i.e. twice the total for random phases. Figure 
[15] illustrates this factor 2 between random and deterministic 
phases. 


V. CHANGE OF MEASUREMENT BASIS 

An intriguing characteristics of multifractal properties is 
their dependence on the basis choice. Indeed, it is known 
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for the intermediate map [|20| that multifractal properties for 
rational 7, which are visible in the momentum basis, disap¬ 
pear in the position basis. It is all the more surprising that a 
recent conjecture |^2l proposes to link the spectral statistics 
(independent of the basis) to the multifractal spectrum (a pri¬ 
ori basis dependent). Apart from its fundamental interest, this 
question is also important for experimental implementations. 
Indeed, it is not always evident to choose the measurement ba¬ 
sis at will in an experiment, and it is thus interesting to assess 
how multifractality is modified when different observables are 
used. The main idea of this section is thus to identify how the 
multifractality spectrum varies when the basis is changed. 

The results of our analysis have shown that in this case the 
multifractality breakdown follows the broad picture of sce¬ 
nario II where the multifractality at small scales is uniformly 
destroyed. However, we have found that this broad picture ad¬ 
mits several variants depending on the presence or absence of 
a characteristic length in the model itself or in the perturba¬ 
tion. In the absence of a perturbation, both the PRBM and the 
Anderson models have no characteristic length, while the in¬ 
termediate map exhibits such a length (see Section II). When 
no characteristic length is present in the unperturbed model, 
like for the Anderson model, we were able to construct two 
kinds of perturbations, which themselves may or may not ex¬ 
hibit an intrinsic characteristic length. 

In the following, we first discuss in Subsection V.A the 
PRBM model and the Anderson model in the case where the 
perturbation does not introduce a characteristic length, show¬ 
ing that in this case the results correspond to scenario II. How¬ 
ever, in Subsection V.B we consider the intermediate map and 
the Anderson model when the perturbation has a characteristic 
length. In both cases, we show through a two-parameter scal¬ 
ing analysis the presence of a perturbation-dependent charac¬ 
teristic length, below which the multifractality is uniformly 
destroyed (following scenario II). However, as we shall see, 
the behavior is different above the characteristic scale. 

A. Absence of a characteristic length 

1. PRBM model 

We first consider the PRBM model defined by (dJ. We con¬ 
struct a generic change of basis through a smooth deformation 
of the identity. The unitary matrix defining the basis change 
is chosen to be 

17(e) = e ieM (54) 

where e is the deformation parameter and M an element of the 
Gaussian Orthogonal Ensemble (GOE) of random matrices. A 
matrix H of the PRBM in the new basis becomes H' given by: 

H' = U(e)HU(e)~ 1 . (55) 

In order to get generic results, we average over a sample of 
matrices M from GOE. 

Fig. [17] displays the curves I)\ (i) of the eigenvectors for 
different e, showing that an appropriate vertical rescaling en¬ 
ables to collapse them (as opposed to the previous horizontal 
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FIG. 17: (Color online) Scaling analysis of the local multifractal di¬ 
mension Di for the PRBM model (j]} for N = 8192 and 6=1 
for different change of basis (155}. Left: raw data. Right: data after 
vertical rescaling of D\ by a factor Di(e). Each curve corresponds 
to a different value of the parameter e in ( 154b . It is chosen to be of 
the form e = 10 0 ' 1 " -3 ; black circles: n*l; red squares: n = 3; 
green diamonds: n = 5; blue up triangles: n = 7; yellow left tri¬ 
angles: n — 9; brown down triangles: n = 11; gray right triangles: 
n = 13; purple pluses: n = 15; cyan crosses: n = 17; magenta 
stars: n = 19. 10 realizations of GOE matrices in d54b are taken in 
order to average over 81920 vectors. 


rescaling performed in Section III). The rescaling here just 
affects the height of the plateau and not the scale at which 
it appears, clearly indicating that scenario II is followed: the 
multifractality disappears uniformly when the perturbation is 
increased. As Fig.[2]shows at e = 0, multifractal dimensions 
are given by the plateau appearing at intermediate scales, and 
therefore the rescaling should be made on these ranges of 
scales. The scaling parameter -Di(e) (corresponding to the 
mean value of the plateau) is displayed in Fig. [18] showing 
that the multifractal dimension smoothly goes to the ergodic 
value for large e. It was checked (data not shown) that the 
same scenario also applies for different values of q. 


2. Anderson model 

Here, we investigate a change of basis for the Anderson 
model. We consider the effects on the multifractal spectrum 
of a rotation U of the critical states of the Anderson model. 
A perturbation of the form (154} would require to calculate the 
exponential of a full matrix M of size L 3 for a 3D system up 
to L = 120. We will use instead the unitary evolution operator 
associated with the so-called quasiperiodic kicked rotor [47]: 

v 2 

Ukr = — + JC(t) cos x ^ 5(t — n) , (56) 

n 

with 1C(t) = K[ 1 + rjcos(uj 2 t) cos(w 3 f)] a quasiperiodic 
kicking amplitude with two frequencies ui 2 = Inty/h and 
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FIG. 18: Scaling parameter D\(e) extracted from the analysis of 
Fig. Q7] It is normalized to correspond to the multifractal dimension 
in the PRBM model {7} with IV = 2 13 , b = 1 in the limit e —> 0. The 
perturbation parameter e quantifies a generic change of basis {55}. 



(a) (b) 



FIG. 20: (Color online) Multifractal dimension D 2 for the Ander¬ 
son model with basis change {57} as a function of (a) t/m associ¬ 
ated to the ID quasiperiodic kicked rotor and (b) t/Txh associated 
with the 3D Anderson model, for different sizes L from L = 40 
to L = 120. The scaling f/Txh put the curves for different system 
sizes on top of each other. However, the ergodic limit D q = 3 is 
reached when t/m ~ 1 associated with the ID dynamics. This 
arises due to the ID character of our rotation. Each symbol cor¬ 
respond to a different size L\ from right to left L = 40 (circles), 
60 (squares), 80 (diamonds), 100 (up triangles), 120 (left triangles). 
Each color correspond to a different value of K = \/2L & Ti 2 /r with 
f = fmaxfT^ 9 /fm^ 9 , n an integer varying from n = 0 to n = 9 
(bottom to top, i.e. black to magenta) and f m ax = L 6 /(21 2 /2Ti 2 ), 
fmin = 50. For each couple of symbol and color, different points 
correspond to different times; t = 10, 12, 16, 21, 27, 35, 46, 59, 77 
and 100 from bottom to top. rxh = N 2 /D and Txh = L 2 / D were 
determined by a numerical determination of the diffusion constant D 
of the ID quasiperiodic kicked rotor. 


FIG. 19: (Color online) Left: Moments P 2 of the Anderson model 
with basis change 157b for different values of t/m, the Thouless time 
associated with the quasiperiodic kicked rotor. W = 16.53 « W c , 
L = 120 and t/m = 0 (black circles), 8.9 10 -11 (red squares), 
1.1 10” 9 (green diamonds), 2.8 10 ' (blue triangles up), 8.7 10 -6 
(purple triangles tilted), 2.3 (brown triangles down). Right: Local 
multifractal dimensions D 2 as a function of t/L for different values 
of t/m, same parameters and color code as Left. In both figures 
each color corresponds to a different couple of t and K ; unperturbed 
Anderson model (black circles); t ss 10, K = 21 (red squares); 
t = 10, K ss 69.5 (green diamonds); t = 21, K = 760.9 (blue 
triangles up); t = 59, K = 2517.7 (purple triangles tilted); t = 100, 
K = 998784 (brown triangles down), rxh = N 2 /D was determined 
by a numerical determination of the diffusion constant D of the ID 
quasiperiodic kicked rotor. 


W 3 = 27rvl3 incommensurate with 27r, K the stochastic- 
ity parameter, t the time, p the momentum conjugate to the 
position x. This ID system is a variant of the famous peri¬ 
odic kicked rotor 1I48L14911 (obtained with 77 = 0), a paradigm 
of quantum chaos known to exhibit the phenomenon of dy¬ 
namical localization, i.e. Anderson localization in momentum 
space. Due to the additional incommensurate frequencies, the 
quasiperiodic kicked rotor performs an Anderson transition 


between a localized phase at small K < K c and a diffusive 
metallic phase at large K > K c {351438} [50]. The evolution 
operator associated with {56} over a unit step of time writes: 


where the value of the effective Planck constant is taken as 
h = 2.89 0. We have used this evolution operator (157} in 
the diffusive metallic phase K K c « 4.7 for // = 0.8 to 
rotate the eigenstates | v k(a)) of the 3D Anderson model {3. 
More precisely, we have considered (i|\E'(a')) = ^(a) as a 
vector of size N = L 3 in the p-space of the quasiperiodic 
kicked rotor and have made it evolve using the evolution op¬ 
erator U over a time t. 

Figure[l9]represents the second moment l\ as a function of 
the box size for different changes of basis, i.e. different values 
of the diffusion constant D of the quasiperiodic kicked rotor 
and different evolution times. In the left panel, P 2 seems to 
scale as a power law of the box size t over the whole range 
accessible. If we study more carefully the local multifractal 
dimensions D 2 (right panel), they show approximate plateaus 
with small variations but no systematic change which could 
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(a) (b) 



FIG. 21: (Color online) Anderson transition for the Anderson 
model with change of basis. Finite-size scaling analysis of = 
lnF^/ In A for fixed A = l/L = 0.1 as a function of disorder 
at various system sizes L £ [20,120]. The results for the standard 
Anderson model are represented in (a) while (b) corresponds to the 
rotated eigenstates with K = 115 and fixed t/Tn ~ 40. The error 
bars are standard deviations. The lines are the fit with a Taylor ex¬ 
pansion of the scaling function T 2 = F(L/£) + L~ v F m (L/£), and 
£ ~ |W — W c \~ u the scaling length. The polynomial fit has ten 
adjustable parameters, v and W c being fixed to their known values 
u = 1.6 and W c = 16.53 (see fl2ll ). In both cases the goodness 
of fit is quite acceptable (larger than 0.T). We find y « 1.8 for the 
unperturbed Anderson model and y ~ 1.5 for the rotated one. This 
compares well with the simulations of fH. 


indicate the presence of a characteristic length. We also see 
that the dimension Z) 2 defined as the average of D 2 over £ 
increases towards the value D 2 = 3 when i/rth —> 1, with 
Tth the Thouless time associated with the quasiperiodic kicked 
rotor, rxh = N 2 /D, N = L 3 and I) the diffusion constant of 
the quasiperiodic kicked rotor. In the case f/rth ~ 1, the 
rotated eigenstates are uniformly distributed over the entire 
sample. 

Note however that different system sizes L lead to a differ¬ 
ent dependence of D 2 as a function of f/rth- When compar¬ 
ing different system sizes, one should consider the Thouless 
time associated with the 3D Anderson model TVh = L 2 /D 
(with D the diffusion constant of the quasiperiodic kicked ro¬ 
tor) instead of tj h associated with the ID quasiperiodic kicked 
rotor. This is done in Fig. [20] Then, the data for Z? 2 are seen 
to collapse onto a single curve, apart from deviations at small 
t/Tj h. This is at the expense of the physical meaning of the 
ergodic limit D q = d = 3 which in this latter case arises 
when f/Txh ~ 10 s . This is due to the fact that we have ro¬ 
tated the eigenstates of the 3D Anderson model using a ID 
diffusive dynamics. This implies a very strong anisotropy in 
3D as the distance between two sites adjacent along the y-axis 
is L in our vector of size N, and L 2 for those adjacent along 
the 2 axis. It is then clear that the ergodic limit is reached 
only when the 2-axis is filled, and this arises when t reaches 
Tth = N 2 /D where N = L 3 is the effective size of the sam¬ 


ple along the 2-axis. 

The remarkable scale invariance observed in Fig. [19] sug¬ 
gests that the rotated Anderson model remains critical at 
W = W c , whatever the amplitude t/Tn, of the perturbation. 
Indeed, away from criticality one expects to observe the emer¬ 
gence of one of the characteristic lengths discussed in IV A 
(the localization length or the correlation length). We have 
checked that the rotated eigenstates perform a localization- 
delocalization transition at the same value of disorder strength 
W c as the unperturbed ones. Following the analysis of 1131 . 
we have considered the quantity T2 = In / In A with fixed 
A = l/L = 0.1 where P% v = exp((lnP 2 )) and (.) stands 
for an average over disorder realizations, see also Sec. IIIBI 
Figure [211 shows the result of a finite-size scaling analysis for 
the standard Anderson model and the rotated one. In both 
cases, we find that our numerical data are compatible with a 
scaling r 2 = F(L /£) + L~ v F- m {L l /'£) where F and Fin- are 
scaling functions and with the localization/correlation length 
£ ~ \W — W c \~ v diverging at W c « 16.53 with the crit¬ 
ical exponent v « 1.6. The irrelevant exponent y controls 
the usual irrelevant corrections (precise values are in the fig¬ 
ure caption of Fig. [21] and quite compatible with the known 
results of j T2ll ). The effect of the rotation is most clearly ob¬ 
served when considering the values of r| = Hm£_j. 00 T 2 (W c ) 
extracted from the finite-size scaling analysis: t' 2 « 1.66 in 
case (a) and t 2 « 2.33 in case (b) of Fig. [2T] The increase of 
T2 with the perturbation strength t/TYh is in good agreement 
with the results of Fig.l20l 


B. Presence of a characteristic length 

1. Intermediate map 

We now consider a change of basis of the random inter¬ 
mediate map <]4]). As for the PRBM model, we construct a 
deformation matrix (l54l > and use it to transform the propaga¬ 
tor U into a new matrix U'. Results are displayed in Fig. [22] 
where the local multifractal dimension D 2 is plotted as a func¬ 
tion of the box size for various rotation parameters e. One 
clearly observes a systematic dependence on l which hints to 
the presence of a characteristic length. However, contrary to 
the results in Sec. Ill, a rescaling of the boxsize is not suf¬ 
ficient to collapse the data. Indeed, a double rescaling (both 
horizontal and vertical) is needed to collapse the data for dif¬ 
ferent e values (see Fig. [22] top). This indicates that there is 
a scaling length £(e) for the box size and a scaling parame¬ 
ter for the multifractal dimension I9 2 (e), both depending on 
e. Above the characteristic length, multifractality is destroyed 
(all multifractal dimensions equal to 1). Below the charac¬ 
teristic length, multifractality survives but the multifractal di¬ 
mensions smoothly go to 1 when the perturbation is increased 
(see e. g. Fig. l22l bottom right). The data are therefore fully 
compatible with a variant of scenario II including the pres¬ 
ence of a characteristic length, below which multifractality is 
uniformly destroyed. Here the characteristic length does not 
come from the perturbation as in Section III A, but originates 
from the intrinsic characteristic length of the intermediate map 
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FIG. 22: (Color online) Top: Local multifractal exponent for the 
random intermediate map, 7 = 1/5, IV = 2 13 , q = 2 under a 
generic change of basis. Left: raw data. Right: rescaled data. The 
points corresponding to l = N/2, N/A and N/8 have been dropped. 
Each curve corresponds to a different value of the parameter t in 
1 1541 . It is chosen to be of the form e = 10°' ln_3 ; red squares: 
n = 4; blue up triangles: n = 8 ; yellow left triangles: n = 9; 
brown down triangles: n = 10; gray right triangles: n = 11; purple 
pluses: n = 12; cyan crosses: n = 13. Bottom: Left: Variation of 
the scaling length £ as a function of t. Right: variation of the second 
scaling parameter D 2 (e) as a function of t. Black full symbols are 
exact data, dashed lines are fitting functions (resp. exponential and 
second order polynomial). 



FIG. 23: (Color online) Discrete Wigner function of an eigenstate of 
the random intermediate map for 7 = 1/3, color/grayness varies 
from blue/dark gray (minimal value) to red/light gray (maximal 
value). The symbols represent lines a\P — ( 12 X = 03 (mod 2N) 
(x = X/2N and p = P/2N). Here N = 128, 0,2 = 1 and ai = 64 
(circles) and 65 (squares). 



I </{(£) 


FIG. 24: (Color online) Left: D 2 as a function of the box size l 
for the random intermediate map with momentum to position basis 
change for 7 = 1/3 and different values of N/ai = 2 s (pentagons), 
N/ai = 2 4 (squares), N/a 1 = 2 3 (circles), N/a\ = 2 2 (triangles), 
N/ai = 2 1 (down-triangles), N/ai = 1 (diamonds), with N = 
2 13 (size of the Wigner function is 2N). Right: Rescaled £>2 as a 
function of the rescaled box size for the same data. 


E (see Section II) which is modified by the perturbation. 

2. From momentum to position basis in the intermediate map 

The intermediate map displays multifractal eigenvectors in 
the momentum representation (Q}. However in the position 
basis eigenvectors are extended. It is natural to ask how mul- 
tifractality is destroyed when one goes from one basis to the 
other. In this subsection we focus on a specific class of ba¬ 
sis changes which interpolate between momentum to position 
basis and are linked to certain types of physical observables. 

The Wigner function of a wavefunction ijj, in our case dis¬ 


crete Wigner function (DWF), is a quasi-probability distri¬ 
bution in phase space and thus provides an adequate testing 
ground to probe the transition from momentum to position. 
We use the DWF as described in m. If the Hilbert space di¬ 
mension is N we define a phase-space grid of 2 N x 2N points. 
Let us label the (2TV) 2 points {X, P), with 1 < X, P < 2 N. 
If the state of the system is given by a density matrix p , then 
the simplest expression for the DWF is 

W(X : P) = Tr[A(X 1 P)p\ (58) 

where A(X, P) are (so-called) point operators defined as the 
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discrete Fourier transforms of the translation operators 

T(ai,o 2 ) = U ai V a2 exp [ina^/N], (59) 

with «i , a2 integers. The translations in ( l59t are defined 
by shifts U and V in position and momentum, such that 
U\X) = \X + 1) and V\P) = \P + 1). In JH it is shown 
that the DWF thus defined complies with all the properties 
expected from a Wigner function. Namely, Yp W(X, P ) is 
the probability \ipx\ 2 associated with the wavefunction in po¬ 
sition representation. Similarly Y xXV (X, P ) is associated 
with intensities in momentum representation. More generally, 
summing W (X, P) along straight lines 

a\P — CL2X = 0,3 (mod 2 N) (60) 

with fixed a±, <12 and varying a3 yields the probability distri¬ 
bution associated with the wavefunction expressed in the basis 
IU3) of eigenvectors of T(ai, o 2 ). In particular for <21 = 1 and 
a.2 = 0 we sum over vertical lines and we get momentum ba¬ 
sis, while for ai = 0 and 02 = 1 we get the position basis. 
Note that since N is a power of 2, whenever 01/02 is also a 
power of 2 the lines will be horizontal for ai > 02 and ver¬ 
tical for ai < 02. We illustrate this in Fig. [23] By changing 
(at, 02) we can go from position to momentum basis. 

In Fig. [23 we show an example of the (absolute value) of 
the DWF for one eigenstate of the random intermediate map 
with 7 = 1/3 as well as two examples of lines (f60l i. 

In Fig. [24] we show the local multifractal dimension D-iU) 
for different values of the slope 02/01 of the lines defined 
in Eq. < f60t with 02 = 1. In this case, the parameter e = 
log 2 (iV/ai) gives the amplitude of the perturbation. The data 
displayed in Fig. [24] show that, as in the preceding case, a 
double rescaling enables to collapse the curves for different e. 
Again, a characteristic length £(e) depending on e separates 


two regimes. Above the scale £(e), multifractality disappears, 
while below this scale it is uniformly and smoothly destroyed 
when e increases (see Fig.l25l>. We conclude that for this more 
physical change of basis the multifractality breaks down fol¬ 
lowing again a variant of the second scenario. 


3. Change of basis with characteristic length scale in the 
Anderson transition 

In the two previous cases considered, a length scale appears 
below which multifractality is smoothly destroyed. In the case 
of the PRBM model however, there is no such feature. We be¬ 
lieve that the difference reflects the fact that the intermediate 
map has an intrinsic length scale 5 = N/b while the PRBM 
model does not. The change of basis rescales this character¬ 
istic length and that is the reason for scaling laws with two 
parameters. 

In this respect, the Anderson model is similar to the PRBM 
model: it has no intrinsic length scale in the critical regime. 
The previous basis change using the quasiperiodic kicked ro¬ 
tor has a characteristic length A = s/Dt associated with dif¬ 
fusion, but since A L for the parameters considered so far, 
this characteristic length does not play any role, as observed 
in Fig.[T9l 

In order to study the effect of such a length scale on the 
basis change, we have to work in a regime where A < L. 
Therefore we consider a rotation with the evolution operator 
of a 3D kicked rotor: 

tt P 2 T .cos x + cos y + cos z v-^ . . 

H 3 kr = y + K - 3 ^-2^(*- n ), (61) 
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FIG. 25: (Color online) Left: scaling length £ as a function of e 
for the random intermediate map with momentum to position basis 
change for 7 = 1/3 and N = 2 13 . Right: second scaling parameter 
Z> 2 (e) as a function of t = log 2 (N/ai). The dashed lines are fitting 
functions (resp. exponential and second order polynomial). Same 
data as in Fig. M 



FIG. 26: (Color online) Double scaling analysis of the local mul¬ 
tifractal dimensions for the Anderson model with isotropic basis 
change. Left: D 2 as a function of the box size £ for different val¬ 
ues of t = t/Trh' from bottom to top, e = 0, 0.000125, 0.00025, 
0.00062, 0.0013, 0.0023, 0.0035, 0.005, 0.0068, 0.0089, 0.011. 
Here Tu, = L 2 / D, L = 80 and D « 0.403 the diffusion constant 
determined numerically for the 3D isotropic periodic Kicked Rotor 
U} with K = 10 and h = 2.89. Right: Rescaled £>2 as a function 
of the rescaled box size for the same values of t. 
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FIG. 27: Left: box size scaling £(e) as a function of e for the An¬ 
derson model with isotropic basis change. Right: scaling parameter 
D 2 (e) as a function of e. Same data as in Fig. 26. 


with K = 10 and K = 2.89. For this choice of parameters, 
the 3D kicked rotor (61) is in the delocalized phase and dis¬ 
plays an isotropic 3D diffusion. We considered the evolution 
over a certain number of periods t of the critical states of the 
Anderson model. We analyzed their multifractal properties by 
considering the behavior of £>2 as a function of the box size £, 
as plotted in Fig. [26] The curves obtained at different pertur¬ 
bation strengths e = t/Tj h depend systematically on £ and e 
and collapse onto a single scaling curve when D 2 is rescaled 
by 772(e) and £ is rescaled by £(e). This strongly suggests a 
picture similar to the scenario obeyed by the intermediate map 
under a basis change: at small scales £ -C £(e), multifractality 
is smoothly and uniformly changed by the perturbation, with 
7?2(e) going from its unperturbed value at e = 0, 772(0) « 1.3 
to 772(e) = 3 when e —> 1 (see Fig.[27l>. 

However, in the present case, £(e) varies similarly to A = 
y/Dt, thus increases as a function of e. In addition, at large 
scales £ £(e), the unperturbed multifractality is recovered. 

In a sense, the multifractal critical states are coarse grained to 
a size A by diffusive evolution, which affects multifractality 
only at small scales (CA. 


VI. CONCLUSION 

In this paper, we have studied the destruction of quantum 
multifractality in the presence of different natural perturba¬ 
tions. The models we considered are representative of several 
classes of systems displaying quantum multifractality. The 
perturbations have been chosen to represent potential experi¬ 
mental constraints in realistic systems. Our numerical and an¬ 
alytical results confirm the conjecture presented in [4]]]. We 
found that multifractality can be destroyed in the presence of 
a perturbation following two scenarios. In scenario I, the per¬ 
turbation introduces a characteristic scale below which multi¬ 
fractality is unchanged, and above which it is completely de¬ 
stroyed. In our case, this describes the smoothing of the sin¬ 
gularity in the intermediate map (Section III) and the Ander¬ 
son model away from the critical point (Section IV A). Sce¬ 


nario II corresponds to a uniform destruction of multifractal¬ 
ity at sufficiently small scales. Depending on the presence of 
a characteristic scale in the system, there can be two variants 
of scenario II. If there is no characteristic scale, multifractal¬ 
ity is the same at all scales and smoothly goes to the ergodic 
value for large perturbation. This is illustrated by the case of 
the intermediate map for a change of slope (Section IV B), 
and by the change of basis for the PRBM model (Section V 
A 1) and the Anderson model (perturbation without charac¬ 
teristic scale. Section V A 2). If the system has a character¬ 
istic scale, scenario II corresponds to the uniform destruction 
of multifractality as before but only below this characteris¬ 
tic scale. This behavior can be revealed by a double scaling 
analysis. Above the characteristic scale, multifractality can be 
completely absent, as in the case of the change of basis in the 
intermediate map (Sections V B 1 and 2), or similar to the un¬ 
perturbed system as in the case of the Anderson model when 
the change of basis has a characteristic length (Section V B 
3). Thus the image presented in lEIfl is confirmed in this more 
detailed analysis, but our results show that subtle variations 
on this broad picture can appear. 

The results presented in this paper also give some insight 
concerning the experimental observation of multifractality in 
various systems. Indeed, experimental setups are unavoidably 
subject to imperfections, which will act as perturbations of 
the ideal model that is implemented, including the measure¬ 
ment in a non optimal basis. The results of Section III show 
that smoothing the singularity of the kicked potential in the 
intermediate map preserves the original multifractality below 
a certain scale, which imposes a minimal resolution to the ex¬ 
perimental measurements. This kind of perturbation can ap¬ 
pear e.g. if the model is implemented with photonic crystals. 
The truncation of a Fourier series to simulate the intermedi¬ 
ate map could be envisioned in a cold atom context, but here 
our results show that a huge number of harmonics are needed, 
and other techniques should be devised to implement the sin¬ 
gular potential. Section IV shows that one can afford an im¬ 
precision on the slope of the potential for finite-size systems, 
but multifractality will be modified, however in a way which 
can be precisely predicted. In an experiment, there will be 
more natural observables corresponding to specific measure¬ 
ment bases. We have investigated the behavior of multifractal 
properties under a change of basis, by interpolating between 
the momentum and position bases, or using a generic change 
of basis built from random unitary matrices. It turns out that 
a modified multifractality is observable for small rotations of 
the basis, but subtle behaviors can emerge depending on the 
presence or absence of a characteristic scale, originating ei¬ 
ther from the model or the perturbation. This is particularly 
striking in the case of the Anderson model, where different 
change of bases could lead to different variants of our second 
scenario. Interestingly, our results confirm that the change of 
basis in this case conserves the criticality of the model. 

Despite many theoretical works in the recent past, direct 
experimental observation of multifractality on quantum wave 
functions has remained elusive up to date. Our results show 
that experimental imperfections will eventually destroy the 
multifractal properties if they exceed a certain level, but that 
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a range of parameters subsists where multifractality could be 
observed. The scenarios for multifractality breakdown con¬ 
firmed in this paper could guide the design of future exper¬ 
iments, which would reliably detect multifractality for small 
imperfection strength and observe the scenarios for larger per¬ 
turbations. 
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